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Abstract 



We advance the previous studies of quantum walks on the line 



C^ ■ with two coins. Such four-state quantum walks driven by a three- 

^ ^, direction shift operator may have nonzero stationary distributions 

(localization) , thus distinguishing themselves from the quantum walks 
on the line in the basic scenario (i.e., driven by a single coin). In this 
^ . work, asymptotic position distributions of the quantum walks are 

^O ! examined. We derive a weak limit for the quantum walks and explicit 

formulas of stationary probability distribution, whose dependencies 
^f^ ■ on the coin parameter and the initial state of quantum walks are 



presented. In particular, it is shown that the weak limit for the 
C3 '_ present quantum walks can be of the form in the basic scenario of 

quantum walks on the line, for certain initial states of the walk and 

certain values of the coin parameter. In the case where localization 

occurs, we show that the stationary probability exponentially decays 

/\ • with the absolute value of a walker's position, independent of the 

parity of time. 



keywords: quantum walks with two coins, three-direction shift operator, 
stationary distribution, localization, weak limit. 

1 Introduction 

Quantum walks (QW) may be described as the natural counterparts of clas- 
sical random walks, governed by the principles of quantum mechanics. Like 
classical random walks, QW are classified into two main types: discrete-time 
QW in 121 13] and continuous-time QW [H [5] . In this paper, we study the 
discrete-time case. Let Z be the set of all integers, the QW on Z, has been 
extensively explored. Among all the findings are the weak limit theorems, 
such as P [HIE]. 



Venegas-Andraca et al. [S] introduced and investigated the QW on Z 
with two entangled coins. For this 4-state QW, its conditional shift op- 
erator embraces three distinct directions. At every time step of the walk, 
depending on the state of the coin, the walker moves forward, or moves 
backward, or stalls at that step. In their numerical simulations of the QW 
[9], the authors found the phenomenon, called "localization", whereby the 
probability distribution of the walker's position is seen to exhibit a persis- 
tent major "spike" (or "peak" ) at the initial position. In what followed, 
Liu and Petulante [10] gave theoretical explanations for the observation of 
"localization" . By the analysis of the spectral properties of the evolution 
operator for the QW, a general formula was derived for the limiting prob- 
ability, from which the authors deduced the limiting value of the height of 
the observed spike at the origin. In this work, we advance the study of this 
model of QW driven by a parameterized coin instead of the Hadamard coin. 
Using the method of the Fourier transform, we show that the normalized 
operator jXt converges weakly to a random variable as t — ;■ oo. The limit 
measure is described by both a ^-function corresponding to localization and 
a density function. As our second main result, we derive explicit formulas of 
stationary probability distribution, whose dependencies on coin parameter 
and the initial state of the QW are presented. A similar treatment of the 
limiting distribution is given by Inui et al. [H] for a 3-state walk governed 
by the Grover coin operator. A closer match to the present model is the 4- 
state QW with a quite different shift operator by Konno and Machida [T2] . 
and Gettrick |T3]. As we will show in this work, the stationary probability 
of the QW is independent of the parity of time t, this is the major difference 
between the present model of QW and the one in [12]. In the case of the 
QW driven by many coins proposed and studied by Brun et al. [11], the 
weak limit measures of the various scaled position operators were obtained 
by Segawa and Konno pLSj. For a one-parameter family of the discrete- 
time QW models in both one- dimensional and two-dimensional lattices, the 
convergence theorems for the moments of the walker's pseudovelocity were 
offered in Ampadu [T6] . Some recent developments with new approaches on 
the asymptotic distributions and localization of discrete-time QW can be 
found in Cantero et al. [T7], Konno and Segawa [18], and Ahlbrecht et al. 

m- 

The rest of this paper is structured as follows: In section 2, we introduce 
notations and definitions for the QW that will be needed in the sequel of 
the paper. Section 3 contains our main results on the limiting distribution 
of the QW. The justifications for the results given in Section 3 are deferred 
to Section 4. A brief conclusion is provided in Section 5. 



2 Formulation of quantum walks with two 
entangled coins on the line 

2.1 Definitions of the quantum walks 

For simplicity, as in |9], the two-coin, or more precisely, two- qubit iTamewoik 
is abbreviated by the name "2cQW" . Analogous to those for the single-coin 
(one-qubit) framework outlined in [201 HI], the coin space of the 2cQW 
framework is the Hilbert space Hec ( "ec" for entangled coin) spanned by the 
orthonormal basis {\j)]j G Be} where Be = {00,01,10,11}. The position 
space is the Hilbert space Tip spanned by the orthonormal basis {\x);x G Z}. 
The "overall" state space of the system is "H = Tip ® "Hec, in terms of which 
a general state of the system may be expressed by the formula: 

X&Z jeBc 

The spatio-temporal progression of the 2cQW is governed by an "evo- 
lution operator" U, which is composed of a coin operator A^c and a shift 
operator S. 

In the 2cQW context, the coin operator is defined as the tensor product 
of two single-qubit operators 

yig(^ ^ j\ Qy Ji.^ 

where A may be any unitary operator acting on a single-coin space. 
The shift operator is given by: 

s = |oo)(oo| (S)J2\^ + i)(^l + |oi)(oi| ® J2 l^^^^l 

i i 

+ |10)(10| ® Y^ \i){i\ + |11)(11| ®Y\i- l){i\. (1) 

i i 

Let / denote the identity operator on Up. Then, in terms of A^c and S, 
the total evolution operator U is given by 

U = S{I(g)A,,). 

Given ipQ E H, where HV'oll = I5 the expression ipt = U^ipo is called the 
wave function for the particle at time t. The corresponding 2cQW with 
initial state ipo is represented by the sequence {V'tjjf • 

Let X denote the position operator on Up, defined by X\x) = x\x) and 

let xjjt = Sa;gz SjeBc ^*(^'-^)l''') ® 1-^) ^^ ^^^ wave function for the particle 
at time t. Then the probability pt{x) of finding the particle at the position 
X at time t is given by the standard formula 



ieBc 



|2 



where | • | indicates the modulus of a complex number. At each instant t, 
the eigenvalues of the operator Xf = W XU^ equate to the possible values 
of the particle's position with corresponding probability Pt{x). 



2.2 Fourier transform formulation of the wave func- 
tion for 2cQW 

As in the single-coin setting treated in |2l [8] , Fourier methods can be applied 
in the 2cQW setting to obtain a useful formulation of the wave function. Let 
\E'j'^(a;) = [il^tix, 1), iptix, 2), ipt^x, 3), il^tix, 4:)]'^ represent the amplitude of the 
wave function, whose four components at position x and time t correspond 
respectively to the coin states 00, 01, 10, and 11. As usual, the superscript T 
denotes the transpose operator. Assuming that the 2cQW is launched from 
the origin, then the initial quantum state of the system is reflected by the 
components of ^^'=(0) = [V^o(0, l),^o(0, 2),^o(0, 3), V^o(0,4)]^ = [«!, «2, «3, «4]' 
where X]j=i 



ai 



1. 



To begin, the spatial Fourier transform of '^'j.'^{x) is defined by 



Akx 



tei 



For instance, under this transformation, the initial amplitude is related to 
its Fourier dual by the formula: 



v^-(fc) = vl/-(0). 



(2) 



In general, the Fourier dual of the overall state space of the 2cQW system 
is the Hilbert space L^(K) ® Tiec, consisting of C^- valued functions: 



0(fc) 



subject to the finiteness condition 






(3) 



nWh'^ 



+ 



^2\\h 



^sIIl^ + II04||l2 < oo. 



Thus, given the initial state '^'^^{k), the Fourier dual of the wave function of 
the 2cQW system is expressed by 

#(fc) = Uec{kf^\k), (4) 

where the total evolution operator Uedk) on L^(K) ® Tiec is given by 



Uecik) 



glfc 
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A. 



(5) 



-ik 



Note that Uec{k) = U{k/2) ® U{k/2), where 

gifc/2 Q 



U{k/2) 







'ik/2 



A. 



Here 



A 



a b 
c d 



(6) 



(7) 



To wrap up this introductory section, we collect some basic facts about the 
eigenvalues and eigenvectors of the operators discussed above. 

Since the matrix A is unitary, we may assume, without loss of generality, 
that its determinant is \A\ = e*^, where ^ is a real constant. Similarly, for the 
unitary matrix U{k/2), we may assume that its eigenvalues are Xi{k) = e*''*^^^ 
and \2{k) = e^^^~^^^^\ where ?7 is a real- valued different iable function of k. 
Let fi(fc) = {vii{k) , Vi2{k))'^ and V2{k) = (t'2i,'y22)"^ denote the correspond- 
ing unit eigenvectors. Since Uec{k) = U{k/2) U{k/2), it follows that the 
eigenvalues of Uedk) are: 



Ai(A;) = [\i{k)f = e*^W where ip{k) 
A^{k) = X,{k)X2ik) = \A\ = e'' 
Asik) = X,{k)X2{k) = \A\ = e'' 
A^{k) = [A2(A;)]2 = e^(2«-^W). 



2r]{k) 



Correspondingly, the unit eigenvectors are: 

Vi{k) = vi{k)^vi{k) 

V2{k) = vi{k)^V2{k) 

Vaik) = v2ik)^vi{k) 

V^ik) =v2{k)®V2{k). 

Finally, in terms of eigenvalues and eigenvectors, the wave function 

\l/f^(fc) may be expanded as follows: 

#(A;) = Ul{k)%'^{k) 

= e'''P^^\Vr{k),%\k))Vi{k) 



+e''\V2{k),mi\k))V2{k) 



+e''\V^{k),%%k))V^{k) 



+e 



it{2e-Lp{k)) 



{V,{k),%\k))V,{k). 



(9) 



3 Limiting distributions for 2cQW 

To prepare for deriving the limiting distribution of 2cQW, we present some 
properties of the evolution operator Uedk). As XDX does in the case of the 
single-coin QW [21j, the expression AZ^A acts as an essential ingredient in 
the formulae for the moments of position probability distribution of 2cQW. 
What follows is the relationship between ADA and XDX, which is easily 
verified based on Eq. ([8]). 



A^DA^ik) =7^^X3 = 
A^DA^ik) = 2\2D\2{k). 

According to ^2]J, we have 



\iD\i + \2DX2 = 0. 



(10) 



:iii 



Eq. (ITT]) implies that f2i'y2r = "^12^2 and f 22^ = "WiiWi- 
In addition, we can show the vahdity of the following identity by a simple 
linear algebra. 



VilVi2 



cXiViiVji - bXiVi2Vi2 



(12) 



For better presentation of our analysis of the limiting distribution of the 
2cQW, we specify the coin operator matrix A as follows, 



^(/3) 



cos (3 sin f3 
sin /3 — cos (3 



(13) 



Here (3 G (0, |). The corresponding evolution operator is restated here for 
convenience. 



g2fc/2 j^Qg ^ ^zk/2 gjj^ ^ 

g-jfc/2 gjj^ ^ _g-«A;/2 pQg ^ 



f/(fc/2) = 

By straightforward calculation, the eigenvalues are 



Xj{k) = ±Jl- cos2 (3 sm\k/2) +icos/3 sin(A;/2), 



so that 



XiDXi = (cos/3cos(fc/2))/(2y'l-cos2/3sin2(A;/2)). 
By Eqs. (ITOl) . (fTTT) and (fT6l) . one can obtain 



cos (3 cos I 



AiL'Ai(A;) = -A4DA4(A;) 



-i/sin^ (3 + cos^ /3 cos^ | 
Eqs. dm), ( TT2I) and ( TT5|) together imply the identities: 

^^ii^^ = tan/3[— J ^ + XiDXi], f2i^^ = -t^iiW2- 

2[Xi + Xi) 



(14) 



(15) 



(16) 



(17) 



:i8) 



According to the methods by Grimmett et al. [5], the moments of the 
position distribution are given as 

/2'r ^ dk 

m\k),D^^r{k))- (19) 

where D = —id/dk denote the position operator in A;-space. 
Using the standard calculations, we arrive at, as t — )■ oo, 

^[(^.An = /"E(^^^)i(^^(^)'^^(^))i'i + o(^")- (20) 

By the method of moments (see [8] and references therein), we may derive 
the following limit theorem. 

Theorem 1. Suppose the 2cQW, controlled by the coin operator A®^, 
is launched from the origin in the initial state ^E'o'^(O) = ai|00) + a2|01) + 
aallO) + a4|ll), where ^,-=i |ajP = 1- For y e [—1, 1], let So{y) denote the 



point mass at the origin and let I(a,b) (y) denote the indicator function of the 
real interval (a, 6). Then, as t — )■ oo, the normalized position distribution 
ft{y) associated with jXt converges, in the sense of a weak limit, to the 
density function 

9 

ft \ X t \ , tan/3/(_cos/3,cos/3)(2/) v;^ ,• ,^^. 

f{y) = coMy) + — ^ — ,. /^ , , = 2^ <^jy'- (21) 

In the above formula, the coefficients cqo, cq, ci and C2 are given by 

Coo = ^ - (sin/3 - l)(|a2p + lagp) + tan^ I3{^j^ - ^^)^Re{ai-al) 

+ (sin 13 — 1) tan /3Re(aia2 + Oiio^ — a2(xl — a^aX) — sin /3Re(a2a3) 
C2 = I — (|a2p + l^aP) — Re(a2a3) + (2tan^/3 + l)Re{a{at) 

+ tan/9Re(a;ia2 + ai«3 — 0:204 — 0304), 
Ci = |aip — |a4p + tan/3Re(Q;ia2 + 0:103 + O2O4 + O3O4), 
Co = I + Re(o203 — 01O4). 

Where 'Re{z) is the real part of a complex number z. 

As an example, we consider the case when /3 = |, oi = 04 = ^, and 
0^2 = 0:3 = 0. Direct calculations shows that the density function in Theorem 
1 becomes 

f{y) = {V2- l)6o{y) + 2y' ^ ^ ' ^ L^ =. (22) 

The component with delta function in Eq. ( 12T|) is an indicator of the 
phenomonon of "localization". In contrast, the density function for weak 
limit measure of the basic model of QW, i.e., 2-state QW, does not comprise 
the delta function. In fact, in the pioneered papers for the treatment of 
weak limit measure of QW, Konno [H [7] and Grimmett et al. [8] obtained a 



limiting density function of the normalized position operator ^ associated 
with 2-state Hadamard QW as follows: 

f(y) = ,, \:'%\ ,, - (23) 

It should be noted that Cqo 7^ is the necessary and sufficient condition 
for localization of the QW on both the quantum coin and the initial state. 
Localization does not always occur in 2cQW. For instance, when /3 = ^, 
ai = 0:2 = 0:3 = — I and 0:4 = |, then cqo = by Theorem 1, and the density 
function of Eg. f l2 1 p just becomes the one given in Eg. (1251) . Interestingly, 
the pattern of function in Eg. (123!) and/or delta function turn out to be 
generic components of the density functions for many existing models of 
QW. Examples of these include, the three-state Grover walks by Inui et 
al. [TT], the QW with multiple coins by Segawa and Konno [13], and the 
four-state QW with a four-direction shift operator by Konno and Machida 
[12] (we point out that this is a closer model to 2cQW considered in the 
present paper, which actually is a four-state QW with a three-direction shift 
operator.) 

Finally, it may be worth noticing that the coefficient Cqo in Eg. fl^Tl) . is the 
sum of the stationary probability limt-^oo Pt{x), i.e., Cqo = "^xez l^i-s-oo Pt{x) 
(see Theorem 2 in [10]). Our analysis of the asymptotic distributions of 
2cQW would not be complete without a consideration of the stationary 
probability of the QW. By using complex integrals as did by Konno and 
Machida [12], we obtain the following explicit formulas for the stationary 
probability. 

Theorem 2. Suppose the 2cQW, controlled by the coin operator A®^, 
is launched from the origin in the initial state ^I/o'^(0) = ai|00) -|- a2|01) + 
ctallO) + a4|ll), where X]i=i I'^iP — 1- Let p{x) = limt^ooPtix). Then the 
limiting probability of finding the walker at \x) is given below. 

(i)p(O) = tan2/3sec2/3(l-sin/3)2(|ai|2 + |a4n + sec^ /3(1 - sin/3)(|a2p + lasH 

+ tan /3 sec^ /3(1 — sin /3)^Re(Q;2a4 + asa^ — ol\OL2 — ai^s) 

- 2tan/3sec/3(l - sin/3)Re(a2a3). (24) 

(ii)p(x) = J+(a,/3)[(l — sin/3)'^sec'^/3]^ for x > 1 and 
(iii)p(a;) = J_(a,/3)[(1 - sin/3)Sec^/3]-^ for x < -1. 

Here J+(a, f3) = tan^ /3[sec^ /3(1 + sin/3)^|ai|^ + |a2p + |«3p + sec^ /3(1 — sin/3)^|a4|^ 

— 2 sec /3(1 + sin (3)Re{aia^ + aio^) + 2Re(aiaI + a2a^) 

— 2sec/3(l — sin/3)Re(Q;2^ + CKsoi)]. (25) 
J-{a,f3) = tan^/3[sec^/3(l-sin/3)Vir + |«2r + las]^ -F sec^ /3(1 -F sin/3)V4|^ 

+ 2 sec /3(1 — sin /3)Re(aia2 + aioi) + 2Re(aia4 + 0:203) 

+ 2sec/3(l + sin/3)Re(Q;204 + O3O4)]. (26) 

Let's revisit the first example where /3 = ^, oi = 04 = ^ and 02 = 
03 = 0, we have p(0) = 3 — 2a/2 by Eg. i^^. This agrees with the findings 



reported in [HI [ID] • Moreover, direct calculations by Theorem 2 show that 
'^xezP(^) = a/2 — 1, which is consistent with the previous result by Theorem 
1. In this case, localization occurs. In the other case when (3 = j, ai = 
^2 = as = — I and 04 = |, it can be deduced that p{x) = for any integer 
X from Theorem 2, thus Cqo = J^x&Pi-^) ~ ^ which is also consistent with 
our previous calculation from Theorem 1. Therefore, localization, in this 
scenario, does not take place. Similar conclusions hold for the 4-state QW 
given in [12]. However, unlike that QW in [12j, the stationary distribution for 
2cQW is independent of the parity of time t, this is one important difference 
that distinguishes one from the other. 

4 Proof of theorems 

This section is devoted to the justification for the limit theorems presented 
in the preceding section. 

Proof of Theorem 1. We begin with the moments of the position distri- 
bution: 

E[{X,/tY] = f^T.(^^y\(VAk),%%k))\'^ + Oit-'). (27) 

By the method of moments (see [8] and references therein), the weak 
limit of Xt/t exists. Let Y be this weak limit. Then we have 

ny<y)=l i2\^v,{kin{k))\'f- m 



i=i 



where h{k,j) = AjDAj{k). 

According to Eqs. (fTOl) and (fT7|) . the probability distribution function in 
Eq. (128!) can be written as 



nY<y) = / m{k),n%m'7r- 

J2arccos(tan/3 , '^ ^ ) ^^ 

/2t ^ ^ ^ dk 

-^ m{kmm?f;^- (29) 

Here H{y) is Heaviside function, which is the cumulative distribution func- 
tion of So{y). 

After taking derivatives of both sides of Eq. (!29l) with respect to y, we 



obtain the density as follows 
dP{Y < y) 



fiy) 



dy 



tan/3 



7r(l — y2)A/r— y^ec^ 

27r 



(yi(fc),V&-(/c))|t2arccos(tan/3^J=^) 



dk 



+ 5oiy) 1^ {mik),%%k))\' + \{Vsik),%%k))n- 



+ 



tan/3 



7r(l — 1/2) a/1 — y"^ sec^/3 



(i4(fc),vi/eifc))i 



fc=2 arccos(— tan /3 



VT^ 



,(30) 



Applying Eqs. ([II]), (Il5]) , (Il6]) and (IIH]) to simplify Eq. (jSO]), one can 
obtain the density function f{y) given in Eq. fl2T]) 



Proof of Theorem 2. By Eqs. (|H]) and f lT5]) . the eigenvalues of Uec{k) can 
be obtained as 

Ai = A^, A2 = A3 = —1, A4 = A2 where Ai and A2 are given by Eq. flT^ . 

By inverse Fourier transformation, the amplitude of the wave function 
of the walker at the position x and the time t is given by 

[7/.t(x,l),V^t(x,2),V^t(x,3),^i(x,4)]^ 



n\x) 



27r 



E 



27r 



-jajfe A t 



K){v,{k),nmvAk). 



(31) 



Using the method of stationary phase as used in pjlj, we may deduce 



that 



n\x) ~ 



/ e-^-\~lfY^{V,{k),n{k))V,{k)^. (32) 



The right hand side of Eq. (!32|) can be converted into complex integrals, 
and the exact values for the complex integrals can be computed. Therefore 
we reach 



nm 



For X > 1, 



n\x) 



'-l\ 



a2-aitanJ'+£^^ 

a3-aitan/3 + 4fi|^ 
. (a2 + a3)tan/3-2aitan2/3 + £^a±k 



1 MZ2 + CI 62-22 + C2 &3-22 + C3 642:2 + 04 



(33) 



„a.'+i 



Zi — Z2 Zi — Z2 Zi — Z2 Z\ — Z2 



For X < — 1, 



,-r,ec/ N / ixt 1 rbiZi+Ci 62^1+02 63^1 + C3 Ml +04, 

^t [x) ~ (-1) ^ttI ' ' ' 



z' 



Z\ — Z2 Zi — Z2 Z\ — Z2 Z\ — Z2 



(34) 



(35) 



Here Zi = — (1 — sin/3)^ sec^/3, Z2 = —{1 + sin/3)^sec^ /5; and 

61 = — («2 + tts) tan /3 + 2a;i tan^ /3 
ci = — (^2 + «3) tan/3 — 2^4 tan^ (3 

62 = («! + «4) tan/3 — 2(0:2 + 03) tan^ /3 + 4ai tan^ /3 
C2 = («! + 04) tan/3 

^3 = (tti + ^4) tan/3 — 2(0:2 + 03) tan^ /3 + 4oi tan^ /3 
C3 = (ai + 0:4) tan/3 

&4 = — (^2 + ^3) tan/3 + (4oi + 204) tan^ /3 — 4(o2 + 03) tan'^ /3 + 801 tan"' /3 
C4 = — (02 + 03) tan/3 + 2oi tan^/3 

The formulas for the stationary probabihty can be derived from Eqs. 
dMD, dSD and ([35D. 

5 Conclusions 

For the quantum walk driven by two coins introduced by [9], we here of- 
fer thoroughly analytic treatments of the asymptotic behaviors of position 
probability distributions. The fact that the walker stalls at each time step 
with a fixed positive probability as indicated by the three-direction shift 
operator, does not imply that the stationary probability must be nonzero 
and localization will surely occur. As shown by the example following The- 
orem 1, localization may not occur for some initial state and some choice 
of the coin parameter; therefore, the density function of the weak limit of 
^ in Theorem 1, in this case, does not contain the delta function as one of 
its components. On the other hand, in the case where localization occurs, 
we show in Theorem 2 that the stationary probability exponentially decays 
with the absolute value of the walker's position, independent of the parity 
of time. Finally, we should point out that Theorem 2 in the present paper 
refines on a previous result (see Theorem 1 in |10]), which also concerns the 
stationary distribution. Indeed, Theorem 2 is a stronger version of Theorem 
1 in [TU]. We wish that Theorem 2 would clarify and correct the comments 
in [To] on the decay rate of the stationary distribution for 2cQW. 
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